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We propose the another, in principle nonperturbative, method of the evaluation of the Wiener 
functional integral for term in the action. All infinite summations in the results are proven to be 
convergent. We find the ’’generalized” Gelfand-Yaglom differential equation implying the functional 
integral in the continuum limit. 


Introduction 

The conventional functional integral calculations rely on the Gaussian method of the integration. This means, that 
physical problems having the second order term in the action are calculable precisely, those problems with higher 
orders terms in the action must be calculated perturbatively. The perturbation methods are applicable in plenty 
of the problems of the contemporary physics, however, more frequently than ever recently we are confronted with 
problems, where conventional perturbation methods are not sufficient. The non-perturbative numerical methods are 
successfully applied to answer the questions in the statistical physics, quantum held theory, if we mentions only the 
well known examples. The numerical methods give the answers ”yes” or ”no” on the quantitative questions, the deal 
of the answer, ’’why” is still missing. 

We are going to discuss the analytical evaluation of the simplest case of the beyond Gaussian functional integral 
calculations, where the action possesses term in the exponent. By the convention, we call the term connected with 
ip'^ the ” coupling constant”. In the conventional perturbation calculations for p'^ theory the results of the functional 
integral calculations are obtained in the form of the power series of the coupling constant (see, i.e., Q). These 
series are asymptotic, divergent, but by sophisticated re-summation procedures one can obtain the reasonable results, 
allowing to take the path integral more seriously. 

In this article we propose the calculation the Wiener functional integral with term of the fourth order in the exponent 
by the another method as in the conventional perturbation approach. In contrast to the conventional perturbation 
theory, we expand into the power series the term linear in the integration variable in the exponent. In such case we 
can profit from the representation of the functional integral by the parabolic cylinder functions. We show, that in 
such case the expansions into the series are uniformly convergent and we find the recurrence relations for the Wiener 
functional integral in the N - dimensional approximation. We find the continuum limit of this finite dimensional 
integral by procedure proposed by Gelfand-Yaglom for continuum limit of the functional integral for the harmonic 
oscillator. 

The article is organized as follows. In part 2 on the example of the one dimensional integral we remember the 
calculation of the integrals with the fourth order term in exponent by the parabolic cylinder functions. In the part 3 
we calculate the N dimensional integral as the approximation of the functional integral. We will prove the uniform 
convergence of the result. In the part 4 we will show on the example of the ’’Independent Value Model” in the sense 
introduced by Klauder Q, that this held theory model we can solve non-perturbatively and we discuss the problem 
of ’’triviality” of this model. In part 5 and in the appendices I and 2 we will study the result of the part 3. In the 
appendix 3 the ” generalized” Gelfand-Yaglom 0 equation determining the continuum limit of the functional integral 
is proven. In the appendix 4 we present the same result as in part 5 for unconditional Wiener measure functional 
integral. 
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Simple integral 


We will define the functional integral as the limit of the finite dimensional integral. This definition has not problems 
with the idea of the integral measure, because for the finite dimensional integrals the integral measure is defined 
correctly. We explain our method of calculation on the example of the one dimensional integral. Calculating the finite 
dimensional integrals we must solve the problem of the calculation of the integral 

+ 00 


/i = 


J dx exp{ —(ax"' + bx'^ + cx)} 


( 1 ) 


where Re a > 0. 

We are interested also to the problems, when the fourth order term is not small, and therefore it is not possible to 
treat it as a perturbation comparing to the rest of the action term. The generally accepted perturbative approach 
rely on the Taylor’s decomposition of the fourth order term with consecutive replacements of the integration and 
summation order: 




(-a)^ 


n—0 


+ 00 

J dx x^"' exp{—+ cx)} 


( 2 ) 


The integrals in above relation can be calculated, but sum is divergent. 

However, R = /i(a, b, c) is an entire function for any complex values of b and c, since there exist all integrals 

+ 00 

b, c) = (-1)’^+’" J dx exp{-(aa;‘‘ + bx'^ + cx)} 


Consequently, the power expansions of Ii = /i(a, 6, c) in c and/or b has an infinite radius of convergence (and in 
particular they are uniformly convergent on any compact set of values of c and/or b). Let us now consider the power 
expansion in c which we shall frequently use: 


oo +°“ 

h = [ dx a;" exp{—(aa;'^ + bx^)} 

n / 

■n.=0 


(3) 


The integrals here appearing can be expressed in terms of the parabolic cylinder function Di,{z), v = —m — 1/2, (see, 
for instance, i)- For n odd, due to symmetry of the integrand the integrals are zero, for n even, n = 2m we have: 


-z’^/4 pco 1 

D_^_ii 2 iz) = — - — / da; exp{—a;^ — za:} 

™ ^ r(m +1/2) 7o 2 ' 




r(m+l/2) J_^ 

Explicitly, for the Eq.Q we have: 

r(i/ 2 ) 


r(m + 1/2) 7o 

+00 1 (_^\zn 


T{m +1/2) 


_ (0 /'z'l f- 

r9„U/4 E D_^_^/2{z) , ^ 


(2a)i/^ 2-^ jji\ 
^ ' 771=0 


4v^ 




•\/2a 


(4) 


(5) 


This sum is convergent for any values of c, b and a positive. 

The convergence of the infinite series m can be shown as follows. For \z\ finite, \z\ < ^/H | arg{—v) \< •7r/2, 

if I oo, the following asymptotic relation is valid [^: 

^ ^ ^ (6) 


(7) 


1 


1 

^ / 1 V 


(ln(-i.)-l)-, 




The m — th term of the sum ([5]) possesses the asymptotic 


— exp 
m! 


—^(Inm — 1) — ^/m z + z^/4 + mln^ 


This means, following Bolzano-Cauchy’s criteria, that the sum ([S]) is not only absolutely, but uniformly convergent 
also for the finite values of the constants of the integral O- 
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The Calculation of the Functional Integral 


In this section we will calculate the functional integral defined as the continuum limit of the finite dimensional 
integral. The continuum limit means the limit N oo, where N is the number of the time slices in the integral. The 
advantage of this method rely in the well defined measure of the finite dimensional integral. Taking the continuum 
limit of the result of the N dimensional integration, we bypass the problem of the continuum integral measure. 

We are going to calculate Wiener functional integral for the continuum action defined by: 


C = 




2 

+ bip{T)‘^ + aip(T)^ 


( 8 ) 


Following the standard procedure, we divide the integration interval to the N equal slices, and we define the N — th 
approximation of the continuum action: 


N 

Cn = ^ A 

i=l 


c/2 



2 

+ + aif'l , 


where A = (3/N. By definition of the integrals we have: 


S = lim Sn ■ 

N—*oo 


We define the N— dimensional integral by the relation Q: 


+ CX5 


N 


Zn = 


n 


dipi 


N 


27rA 


exp<^-^A 


/ri / ‘fii — l 1 I I, 2 I 4 

c/2 ( -^-) +bip^+ a^p^ 


The Wiener unconditional measure functional integral is defined by formal limit: 

Z = lim Z]si . 

N —»-oo 

We use the formal notation for the continuum Wiener functional integral: 

P 


2^ = y \Vlp{x)] exp 


(9) 


The most important problem is to calculate the dimensional integral ®. We rewrite the action in the convenient 
form for the consecutive integrations: 


4 c , &A^ 2 c 


Cn = Aaipf + ^-(l +- )‘pi - -r^i^2 

Ac A 


A 4 C &A^, 2 C 

+ Aaifi + —{l + ——)Lpi - 


( 10 ) 


A 4 c 6 A^^ 2 

+ AaipN + - )Cn 

A c 

In the integration over variable (pi one find the terms resulting from the precedi^ integration over variable Pi-i-, 
as one create the terms for integration over variable pi+i- Exploating the formula [6|: 


y ^ exp(-pa:^ - qx) dx = T{a){2p) “A exp . 


( 11 ) 


for integration over the variable (pi first, we find: 
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expj-Aatpf - + ■ (12) 

By Taylor’s expansions of the term 

exp 

taking into account that terms of the odd powers in 931 disappears due to the symmetry of integral, we simplify the 
integral by substitution: 

tfl = to. 



We reads: 


+ 00 

Z, = ^ f f 1 (u;)'=-i/^exp^ -Aatu^ 


ki —0 

Using formula m we find: 


{2ki)\ 


2-n-A 


c bA"^ 

A<‘ + —>“ 


( 13 ) 




fci=0 


(2fci)! I / 2nA 




fel+1/2 


where 


z = 


c(l + 6 A^/c) 
v'2aA3 


Following definition of z, we can use the factor 


to simplify further this relation to the form: 


A 


(1 + 6 A Vc) 


fci+1/2 


exp(zV4) D_k,-i/2{z) , (14) 


= 


{^ 2 )^ ^ f c/A 


—( 7 - , ) r(fci + 1/2) exp (z^/4) D_i. i/ 2 {z) , 

V2^(l + 6AVc) {2k,)\ \{l + hA^lc)) ^ 

or, equivalently: 


(15) 


= 


E 


( = V 

V A(l+bAUe )) 


(¥’ 2 )"'=^ r(fci + 1/2 )I?_,,_i/2(z) , 


\/27r(l + 6A2/c) 
where the notation was used: 

V_k,-i,2{z) = z'^^+i/^eX D_u,.y2{z) ■ 


(16) 


(17) 


From above result, the term ((^ 2 )^^^ will be involved into the second integration over the variable tf 2 - This will results 
to the additional fci dependent contribution to the result of the first integration and to the introduction of the link 
between both steps of integration. Explicitly, for the second integration, we have: 


+ 00 


= 


dtp2 


2ttA 


{‘P 2 ) 


2 ki 


exp 


r hA‘^ r 

-Aa^t - 7-(l +-)7’2 + 7-<^2(^3 

Ac A 


(18) 
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By Taylor’s expansion of the term of the action linear in variable ip 2 -, by the interchange of the order of the integration 
and summation, taking into account that the odd powers of tp 2 variable disappear due to symmetry of integral, and 
by variable substitution = <^ we find: 


4°^ = E 




+ CXD 


(2fc2)! J 

k2=0 ^ 0 

By the same evaluation as for Z\ we have: 


' H'=^+'=i-i/2exp|-Aaa;2 - -^(i + ^ . 

2ttA I Ac 


(19) 


'jloc _ 

^2 — 


(c(l + 6AVc)/A) ~ 


k2 


\/27r(l + &A7c) 

Taking both integration steeps together we have: 


f _ c 


^2 = 


E 




k2 


( 20 ) 


\\/2’r(l + 6A2/c)^ 
where the new symbol was introduced: 


(27)1(27)! 


r(A:i + l/2)I?_fc^_i/2(z) r(/ci + 7 + l/2)T’-/ci-fc2-i/2(^) ; 

( 21 ) 


(1 + 6A2/c) 

We can repeats this procedure for the integration variables - For the integration over variable pN, 

there one have no linear term in exponent, therefore we don’t expand anything and this last integration don’t add 
the summation to the final formula. We have: 


-■N 


+ 00 



— CXD 


dpN 





exp 


— Aa(f 


4 

N 



( 22 ) 


and the result is: 


(c(l/2 + &A2/c)/A)"^-^ 


'^Loc _ 

— 


r(fcAr-i + l/2)'D_k^_^_i/2izN) ■ 


^2'k{1I2 + bA'^/c) 

Repeating by recurrence this procedure, we obtain finally the non-perturbative, exact result: 


2 /„M-1/2 


oo N 


Zat = [27r(l + 6AVc)] [27r(l/2 + feA^/c)] 


E n 

ki,--- ,kN — i—0 i—1 


(?.) 


2ki 


(23) 


(2fc, 


—r(7-i + 7 + i/2)i2_fe._^_fe._i/2 (zi) 


(24) 

where 7 = = 0, and 7 = 7 = ’ ‘ ‘ = ^n -2 = Cj = Q also zi = 2:2 = • • • = zj^-i = z. The modifications appears 

due to the last integration step and we have 


c(l/2 + 6A^/c) 

ZN = - , ^ „ - , ?Af-l = 

y/2aA^ 


1 


1 


(l + 5A2/c) V (l/2 + 6A2/c) 


For the conditional measure functional integral, the variable pN will be fixed, in time-slicing procedure the last integral 
over variable pN is absent. This means, that we’ll have the same values for all 7 as well as for Zi, and upper limit for 
product formula in (l24l) should be iV — 1. 

In the noninteracting harmonic oscillator limit, {b positive, a = 0), the Eg. ((Mil is reduced to the TV-dimensional 
approximation to the unconditional measure Wiener functional integral for the harmonic oscillator. For a 7 ^ 0 we can 
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derive from (l24)l by Gelfand-Yaglom procedure the differential equation of the second order. The dependence on the 
coupling constant is in the part characterized by multiple summations over indices ki. This part is now the object of 
our discussion. 

We are going to show that relation (IMl) for the dimensional integral is uniformly convergent for each of the 
summation index. For the summation over one index ki we have (we omit the index i in ^i and Zi)\ 


E 


(0 


2ki 


{2ki)\ 


r{ki_i + ki + (^) 


[r(A:i + fci+i + (-z)] • 


(25) 


For ki ^ oo we apply for the D functions the asymptotic relation ([HI). We have for the ki dependent part of (1^51) the 
relation: 


O-ki = 


(^)2 tY[ki + ki+i + l/2)T(ki-i + ki + 1/2) 


exp {z^j2) 


(26) 


X 2;fci+fei+i + l/2 gxp 


X + gxp 


(fc^ + ki+i) - 1) - \/ki + ki+i z 


ki-i + k 


(In (fci_i + ki) - 1) - -yfelT+l)/ 0 


We convert the above relation by help of the identity: 

r(p + fc + 1/2) = exp [lnr(p + fc + 1/2)] 

to the form convenient for the proof of the convergence. 

For the logarithm of the relation (l26l) we have: 

Inofe; = lnr(A:i + fci+i + 1/2) + lnr(fci_i -\- ki + 1/2) — lnF(2A;i + 1) 

- (fci_i + ki) - 1) - EiE±i(ln [ki + h+i) - 1) + 2k^ In(^z) + (fci_i + ki+i + 1) Inz 

~ i^/ki + fci+i + -s/ki + ki-i)z + z^ 12 + o(fcj ^ ) 

For the gamma functions with argument u ^ oo we can use the asymptotic relation : 

InF(u) = (u — 1/2) ln(u) — u + l/21n(27r) + 0(u~^) 

For ki-i and fci+i fixed and finite and if < ki, we find the following asymptotic for ak^- 


(27) 


(28) 


Inofe, = -{ki- 


ki-i + ki+i — 1 


) Infcj + 2fcj(ln(^z) — ln2 + 1/2) — 2-\/^ z + (fci_i + ki+i + 1) Inz 


+ 1/2 InTT + z^/2 + o{k^ ^^^) 

The leading term of the above relation is 


(29) 


-fciln(fci) , 


and Gki is going to zero in the asymptotic region as 


kffi^' 

ki\ exp {^/h) 

where a and (3 are finite numbers. The asymptotic feature of is sufficient for proof of the uniform convergence of 
the series (051) and therefore for the finiteness of the sum (1^5)) . 

By the same method we can prove the uniform convergence of the summation over two, three, ..., — 1 summation 

indices ki in the equation (12411 . This important conclusion indicate, that the N dimensional approach to the the 
continuum functional integral can be summed up. 
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Independent Value Model 

The another example of the use of the calculation of the functional integral by help of the parabolic cylinder function 
is the problem of the independent value modelQ- The independent value models (ivm) are covariant models without 
gradient terms. The formal functional integral for generating functional of such models reads: 


=J\f J Vcj) exp|y J{x)(j){x) — -w?(j){xY — 5 (?!)(a;)^]| , (30) 

where <i){x) is a real field, J{x) is the source term and g is the interaction constant. The normalization constant M 
guarantees that = 1. It should be noted, that physically the model is trivial (the field equation reduces to 

4>{x) = 0), corresponding to independent fluctuations for different x’s. However, the functional integrals of this type 
have a well-defined mathematical meaning, and its investigations can serve as a tool for a better understanding of 
various approximative schemes of calculations. 

The generating functional of the ivm can be calculated as the continuum limit of the discretized finite 

dimensional integral approximation to the continuum functional integral 

N N ( N \ 

n = M) / d(j)k exp i '^{iJk(j)k^ - l/2ml(j)l/S - go^feA) \ , (31) 

k=l k=i U=i J 

where A is the volume of the discretized space cell, and Aq is the normalization constant. The values toq and go are 
supposed to be both positive and A dependent. In this model the A dependence is chosen so that in the continuum 
limit A —> 0 the physical observables in the model take reasonable values. Only when this dependence is fixed the 
model is defined. The lnZ’[J]“’" is the generating functional for the connected correlation functions and all even 
order connected correlations functions are nonnegative, as follows from the standard canonical form of the generating 
functional for the ivm model. Performing a discretized calculation Lebowitz Q proved that the connected correlation 
functions of the fourth order are non-positive, when a free measure of the functional integral was used. The two 
inequalities for the connected four point correlation function are simultaneously valid only if such function vanish. 
This means, that the discretized approach to the functional integral leads to the gaussian generating functional of a 
free, noninteracting model, with no dependence on the coupling constant in the continuum limit. Such behavior is 
known as ” triviality”. 

We show that we can obtain a nontrivial result for the generating functional of ivm model non-perturbatively by 
taking a particular N —> oo limit in (|3ip - the continuum limit of the discretized approximation. In evaluations of 
finite dimensional integrals we follow the parabolic cylinder functions representation of the integrals with x'^ in the 
exponent of the integrand: 


+ 00 


x^'^dx exp{—(aa;'^ -I- bx^)} 


= + 1/2) exp (^) (^^) . (32) 

where Re a > 0, and D,j{z) is the parabolic cylinder function of index v and argument z. 

Applying (1321) to the finite dimensional approximation of the functional integral (|31D . we find for the k — th factor in 
m the result: 


OO /_ \2 2 

zr™[j] = (Afo)'/^ E exp (^)4 D_._i/2 {z) , (33) 

i=0 *■ 

where pk = (J^A)/(2mQ) and z = {rriQ^/S)/y/8go- The sum in Eg. (1551) is uniformly convergent with respect to z 
on any compact domain as follows from the asymptotic expression for parabolic cylinder functions when the index 
\v\ oo. 

The definition of the discretized approach is accomplished when the dependence of ’’bare” parameters in eD 
on A (or equivalently on N) is fixed in such way that continuum limit N —> oo is well-defined. The requirement 





limA^o(5o/'^) 7^ 0 (i.e., z ^ const in the continuum limit) is the necessary condition for a finite continuum limit of 
the discretized ivm model. Using (1551) the exact formula for the generating functional (1511) reads 


V] 



- PkZ 


D-3/2{z) 

D-ii2iz) 


1 2 2 


D-5/2{z) _ 1 \ 

D- 1 / 2 {z) " J J 


(34) 


where we have taken into account the proper normalization = 1. In the continuum limit ^ A is replaced by 

f dx". We see immediately that the continuum limit survive only terms with the number of the summations equal 
to the power of A. This restricts the contributions to terms up to the first order in pk- In the limit N oo (i.e., 
A ^ 0) we obtain the nontrivial contribution to the generating functional (I30p when mo and ~ pqN are kept 
fixed: 





d'^x 


JHx) ] 

2ml j 


k(z) = I — 


^^-5/2(z) 

2D_,/2(z) ’ 


(35) 


The result ((35l) represents a free theory result, due to the quadratical dependence in the exponent. However, 

()35p possesses a nontrivial dependence on interaction term which survives in the continuum limit and is hidden in 
z ~ yqA/go) dependence of the factor k{z) multiplying the integral in the exponent. Introducing a new parameter 
w?{z,m^) = all non-trivial dependence on z and is included into mass renormalization. This is a new 

feature, because within perturbative approach was found no dependence on the interaction constant in the continuum 
limit. 


Summation over ki 


The result of the Wdimensional integration (1241) is the exact one, calculated without any approximation. The 
multiple summations suppress this advantage somewhat, therefore we shall attempt to provide the ki summations 
in the formula (IMl) . The details of this calculations are in the Appendices I and 2, here we sketch the method and 
results only. 

To obtain the formula for subsequent evaluation of the functional integral, we must to solve the problem of sum¬ 
mations over indexes ki in Eq. (1241) . For each summation index ki we meet the sum of the series with product of two 
T) functions. We are dealing with the sum of the series: 


2ki 


o-ki — 0! + ki + I/2)I?_^,._j_j,._i/2 {z)T{ki + fci+1 -I- I/2)I?_/j._fe._|_^_i/2 {z) 


(36) 


Of course, parabolic cylinder functions belong to the representation of the upper-triangle matrices, therefore we 
believe to the simplification of the product of two such functions. What follows, we propose the our approach to 
provide the summation over indices ki. At first, we represent one of D functions by Poincare - type expansion [^, 
valid for real index and positive argument of the function, if the inequality a < z is valid: 


V4 ^a+1/2 D_^_,/^{Z) = V_^_,/^{Z) = (-1)^' 

1=0 


(g -I- I/2)2j 

j\ (2z^y 


+ £j-(a,z) 


(37) 


where ej{a, z) is the remainder of the Poincare - type expansion of the V function. 

Inserting this relation to the sum, we divide the sum over ki into two parts. One, over finite where Poicare - 
type expansion is correct and the second, the remainder over infinite ki indexes, but small comparing to first part due 
to uniform convergence: 


(—I)-^ 


^2 2 \j (' 2 ^.')|^(^i-i + + l/2)r(fcj -I- fci+i -I- 2j + I/2)T>_i,._j_fe._i/2 (z) -I- 

j=o d- \ z ) (^ i). 

+ y]] 7^Tp^r(fci_i-I-fci-I-I/2)r(fci-I-fci+i-I-I/2)T>_fc._j_/j._i/2 (^) ej(«i 2 ) + 

fei=0 ' 


( 38 ) 


+ E jfEr(fc,_i + fc, + i/ 2 )r(fc.+ ki+i + I/2)(D_fe._j_fc._i/2 (z)I?_fc._fe._^^_i/2 (2), 

U. —M^_L1 ^ 


2ki 


ki —iVo + 1 
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where Nq < z. In the first term of the above relation we add the terms expanding the summation over index ki to 
infinity and in the as ymp totic region of ki > Nq we expand the function V by double asymptotic properties expansions 
proposed by Temme [ll|: 


exp i-Az'^) 

where the following quantities were introduced: 


E 

.A;=0 


fkiX) 


y2k 


y2n 


Rn{a,z) 


(39) 


A = , Wo = - [Vl + 4A — 1] , ^ i Wq + wq — A — A In Wo + A In A . 


The functions /fc(A) are calculated in [HI- For sum of the series (IM)l we find: 


J 


E E + l/2)r(fcj + fcj+l + (z) + 




^ — (0 
ki=0 


2ki 


(40) 


No 

E I'or. M + ^* + l/2)r(A:i + {z) €j{a,z) + 


ki=0 


i2hy. 


E 


iO 


2ki 


ki — 


ki —iVo + 1 


(2A:,)! 


^ .^exp (—ylz^) 
(l + 4A)i/4 

■^Efe(A)' 

.k^O 

^ ^exp(-y42;2) 

1/2 (^) (l + 4A)l/4 

1 

■^Rn{a,z 


J 


\2ki 


(-1)^ V 

'! {2z'^)^ 

i=0 ^ ’ fci=Afo + l 

The sum over index ki in the first term of the above relation will be provided by help of the formula 


“ E -I (2 ‘ 2\3 E 7^T^r(A:i_i + fcj + l/2)r(fcj + fcj+i + 2j + 1/2 )T>_j.._j_/;._i/ 2 (z) ( 41 ) 

o—n J' \ ^ ) 1..— AT _LI A «/' 


e"'^^E fc(^) = (x - t) . 


(42) 


Method of this sum is described in Appendix 1, we are going to evaluate the three terms of remainder now. First of 
all, we must to develop the criteria for the remainder. Schematically, we must evaluate the multiple sums: 


^ ' Ofcl,fc2®fe2,fe3 ' ' ■ ^kM-ltkffO'kN ,kN + l 

ki ,---,kAr—0 


(43) 


If each individual sum exist and it is finite, we divide the sum over index ki to finite, principal sum and the remainder, 
which can be done as small as possible, comparing to the principal part of the sum. Let us define the principal sum 
object, where sum over m indexes ki was not provided, but for indexes km+i,---,kN the finite sum was done: 


^m,N—m — ^ ^ G,ki,k2^k2,ko ‘ ‘ ‘ ^km-l,km^kmi 

kl j-'-ikm—O 


(44) 


where 


No 

^ ^ ®fcm.fem + l®fem + l.fcm + 2 ■ ■ ' + l (45) 

km + l','",kN—0 

The next sum over km we divide into the finite, principal sum and the remainder, the infinite sum, giving minor 
contribution. We have the identity: 


Alyn.M—m — ^m — l,M—m+l T — — m 


(46) 
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where 




CXD 

km — 


Let e = max (ei, £ 2 ,.. •, cat), then we find: 


Sat^o < So ,at + 


So,A/-l £ + ••• + 


So.Af—i £* + ••• + £ 


N 


If Eo,fc are the members of non-decreasing series, the following inequality is valid: 

So,Af < Eat.o < So,A/(l + 


(47) 


(48) 


In the limit N —> 00 we obtain So,oo = Soo,o, if £ ~ N~^~^ where e > 0. Let us check if the terms in the remainder 
corresponds to this demand. 

For Poincare - type expansion the upper bound of remainder was calculated by Olver [^. We use the improved 
upper bound evaluated by Temme [T^. The upper bound for remainder in definition ([37)l can be read: 


ej-(a,z) I < 


z2 - 2a (J - 1)! (2z2y i^ 2(J/2,1/2;J/2 + 1-,1 ^^)exp 


4 .A 

iF2(1/2,1/2;3/2;1-^)) (49) 


— 2a 


This relation is valid for 2-^/a < z. The following quantity was introduced: 


<5 = 


3 

T 16 


2a / a 

+ ^ 




2z^J (z2 - 2a)2 


The sum in remainder: 

NIr(^i-i + ki + l/2)r{ki + ki+i + l/2)I?_fc._j_fc._i/2 {z) ej{ki + ki+i,z) 

j.._n 


(50) 


ki=0 

is bounded by the relation: 

M 


{J - 1)! {2z^y {2h)\ 


^0 /^\2ki 

+ l/2)r(fci -I- ki+i + 2J + 1/2)1? -fei_i-fci-l/2 (^) 


where 


M. = max 


-1 F2{J/2, 1/2; J/2 + 1; 1 - J) exp 


46 , 11 , o? ^ 

iF2(1/2,1/2;3/2;1--) 


— 2a 


(51) 


(52) 


and a= ki + ki = 1,2,... ,Nq . We have the freedom to choose the parameter J, therefore the upper bound on 
this contribution to the remainder can be done as small as necessary. 

For the remainder in the double asymptotic property expansion fllj | we find the definition 0: 


Rn{a,z) = 


,,2a+l/2 


r(2a -f 1/2) jg 

It can be shown by some algebra, that upper bound on remainder is: 

Rn{a,z) < Mn , 


/„(s)^ 

vs 


(53) 


where M„ is upper bound to the function /n(s). 

We are now going to estimate the upper bound of the following contribution to the remainder in (j41ll : 


E 


(0 


2ki 


exp {—Az^) 


2^ r 


ki —A^o + l 


(2fc, 


■yyF(A:i_i + ki + l/2)r(fci -I- ki+i + l/2)I?_fe._^_fe._i/2 (z) ^ 4 ;y)i /4 




Rnia, z) 


(54) 


a—ki-\-kij^\ 
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In the asymptotic region of ki we use the Stirling formula for the gamma functions: 

Inr(u) = (it — 1/2) In u — It + 1/2 In (27r) + 0{u~^) 

We find: 

fci+i+fei_i-i/2_ f fci+i(^i+i + 1/2) + ki-i{ki-i + 1/2) 


(55) 


(2k )l ^(^*“1 + + l/2)r(A:i + ki+i + 1/2) = 


ih) 


exp 


ki 


(56) 

Moreover, in the asymptotic region of ki we replace the function (z) by the principal part of the double 

asymptotic expansion: 


V 


-ki-i-ki- 1/2 


iz) = 


exp {—Az^) 


(l + 4A)i/4 


Taking into account that 




2 V * ' l + \/l + 4A 
and that for A > 0 the following inequalities are valid: 

2 . 2 


+ Az^ In 


1 + Vl + 4A 


(57) 


(58) 


1 - 


1 + Vl + 4A 
we find the upper bound to remainder (1541) : 


< 


“ ^ ’ 1 + VTTIx ~ Vx ’ 


^2n / ^ 

ki—NQ-\-l 


2ki 


(h) 


fei+i+fei_i-i/2 /fct+i(fci+i + 1/2) + fct_i(A:t_i + 1/2)^_^ ki-i-i + ki-i 


exp 


exp [ ki -\- 


h 


1 /^ y,‘2ki-\-ki^\-\-ki — \ 


ki^\-\-ki 


ki_\-\-ki 


{ki) 


ki + - 


1 + 


1 + 


ki-1 


(59) 


For ki big, we use the identity: 


ki±i+ki 


1 + 


fci±i 


1 - 


4i±l 

ki 


ki±i+ki 


/ ki^i , 

exp(-^) 


(60) 


Finally, we have: 




2ki 


_ y" (k) 

j2n-fci+i-fci_i-l/2 I 1 2 / ^ 

ki=No + l ^ 


>=i +i+>=i-i 1 /fci+i(A:i+i + 1/2) + fci_i(fci_i + 1/2) . 

exp ( -- ) (61) 


The function 


exp 


ki+i{ki+i + 1 / 2 ) + ki-i{ki-i + 1 / 2 ) 


is decreasing for ki > Nq, therefore following mean value lemma we read: 

2k 

^271 —— 


^^-73175 exp (p.) E ^(yj 

ki^No + 1 ^ ^ 




where pk = {ki+i{ki+i + 1/2) + ki-i{ki-i + l/2))/2, fc e (iVo,oo). The series 

Nq + 1 


k\ \ A 


(62) 


(63) 


(64) 


































is the remainder of the series uniformly convergent to the function 

fix) = ixd,)^ e^ 


12 


yi 2 

where x = and m > + ki-i — l)/2. By choosing the parameter n, which corresponds to number of the 

decomposition terms taken into account of double asymptotic expansion of the function V, this contribution to the 
remainder can be done as small as we need. 

Now we are going to estimate the last part of the remainder, corresponding to the difference of the series: 


E 

ki —A/^q + 1 

J 


( 2 ! k + l/2)r(fci + fcj+i + l/2)'D_k-_j^-ki-i/2 (z) 


E 

.k=0 


fkiX) 


y2k 


-E 


(-1V 


E 


(0 


2ki 


7! (2z2)i ^ ( 2 /cN! 

7=0 ^ ^ ki=No + l ^ ' 


r(fci-i + fci + i/2)r(fci + fci+i + 2j + (z) 


(65) 


Following the identity (1561) we can recognize in the sums over ki the remainders of the uniformly convergent series. 
These remainders are small, if A = A:? jz is the large quantity. Due to the uniform convergence, we can change the 
order of the summations and we read: 


00 


E 

ki — 


(2fc,)! 


r(fc,_ 


1 


ki + l/2)r(fci + ki-\.i + l/ 2 )I?_fe._j_fe^_i /2 {z) 


exp i—Az'^) 
(l + 4A)i/4 


'n—1 


E 


fkjX) 

z2k 


i^yih + h+^ + 1/2)^^ 

j=0 


jl (2^2)7 


( 66 ) 


The double asymptotic expansion reduces to Poincare - type expansion [T^ when a is fixed after expanding the 
quantities in expansion that depend on A = a/z^ for small values of this parameter. In the difference 


exp {—Azy 
(l + 4A)i/4 


E 


. fe =0 


fkjX) 

z2k 



7=0 


i—l)yki + ki+i + 1 / 2 ) 


27 


j! (2 z2)7 


(67) 


all terms where fc, j < min (n, fl) cancel one another and the rest terms are proportional or smaller then A,j) ^ 

As in the case of previous contributions to the remainder, this means, that this part of the remainder can be done as 
small as we need. 

We shown, that we can evaluate N— dimensional integral (l24)l by method of Poincare decomposition of one of 
the parabolic cylinder functions. We are able to accomplish the sum over the summation index of the Taylor’s 
decomposition of the linear exponential function. We shown, that in the continuum limit the N— dimensional 
integral evaluated by this method have the same value as the corresponding continuum functional integral. 

The result of the recurrence summation procedure of the principal contribution is (see Appendix 1): 



f ^ 

1-'/“ J 

Zn = i 

n [2(1 + 6AVc)u;, 

1 E 



) 

Symbol iN)2j i is defined by the recurrence relation: 



/ \ 

2A 

(«+ 1 ) 2 m.P = E 1 


E 

A=0 

\ / 

i—max 

[ 0 , p-2(/i-A)] 


(A = i)2.,p = -It 

UJq 


(- 1 )^ 
y. (2^2)^ 


(TV) 


2fL, 0 


UJa — lUJct 


(«)2A.* 


yA-A+i 


where 


iOi 





( 68 ) 


(69) 
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A 


UJQ 


1 


2(1 


bA^/c) 

1 


(l/2 + 6A2/c) 


(70) 


-Zjv in the relation (I68|) is the N — th approximation of the functional integral. The continuum limit oi Zjy we 
calculate by the procedure proposed by Gelfand- Yaglom for the calculation of the functional integral for the harmonic 
oscillator [^. The method and the solutions are discussed in the Appendix 3, here we present the result only. 

The functional integral in the continuum limit is defined by the relation 


lim 


AT—>oo 


2n 


1 

7W) 


where the function F(P) is the solution of the differential equation taken in the point /?: 


d 


d 


d 


2b 


92 


—A(r) + 4-Fir) - In 5(r) = F(r) - - 2— In Sir) - 4 - In Sir) 


9r2 


dr 


dr 


dr^ 


d 


dr 


(71) 


The initial conditions are: 


The function S'(t) is defined as: 


F(0) = 1, 
|:F(0) = 0, 


(72) 


(73) 


The equation GB) can be simplified by substitution 


Fir) 


y{T) 

{Sir))^ 


For j/(t) we find the equation: 


^,y{r) = yir) - 


(74) 


In the present time the calculation of the function S'(t) is not finished, therefore we present the calculation of the 
lowest order term in the power expansion in the coupling constant in the Appendix 2. For the first two terms of such 
expansion in the continuum we find: 


lim y 

A^O 2^ 

fi—0 


(- 1 )^ 

(2z2)/^ 


iN) 


2 ( 1 , 0 


= 1 - 


( 1 / 2 ), 


(. 2^3 


{tanh(T 7 ) + ry [3tanh^(T7) — l] } 


where 



Conclusions 

In the present article we presented the new method of the calculation of the Wiener unconditional measure functional 
integral for the action with the fourth order term. This method can be extended to the case with n = 3,4,... 
term in the action, but we did not discuss these possibilities here. The main results are the analytic formula for the N 
dimensional integral and the generalized Gelfand-Yaglom equation implying the functional integral in the continuum 
limit. 

Acknowledgements. This work was supported by VEGA projects No. 2/3106/2003 (J.B.) and No. 1/025/03 
(P.P). The authors acknowledge J. Polonyi and M. Znojil for the discussions and interest about this work. 
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Appendix 1 


In this appendix we propose and on pedagogical level to evaluate the ki summations by the reccurence method. 
Let we start with summation over the index fcjv-i of the Eg. (12411 . The sum that must be done is: 


E 

kN-i—O 


/ ^2 \ f^N-1 


r(fc^_i +1/2)1? — fciV-i —1/2 (^o) 


\/27r(l + blS?/c)ijjQ 

Let us remind the dependence of the function 2? on the parabolic cylinder function D\ 

I?_„_l/2(z) = D_^_y2{z) 

and, also the definitions of the variables loq, z and zq\ 

c(l + bA^/c) 


z = 


zo = 


e = 


Wo = 


^/2aA^ 

c(l/2 + bA'^/c) 
Y^2aA3 

1 


1 + bA^ jc 
1/2 + 6A7c 
1 + bA'^fc 


(75) 


The sum (j75ll is uniformly convergent, therefore for an arbitrary small positive number e exist the number Nq such 
that the following inequality is true: 


OO 

I E I < ( 76 ) 

/CAT-l—AT q + I 


When we replace the true sum by the truncated one, we introduce into our calculation error up to the order N~^~^. 
All recurrence procedure consists of the N steps, therefore we require that in the continuum limit the influence of the 
reminders of the sums disappears. For the truncated sum we demand that 

No < z ^ 

Therefore we are permitted to use the asymptotic Poincare type expansion of the parabolic cylinder function, which 
for 2? means: 


^-fcn-l/2(^0) 


^fcn + 1/2 
^0 


o/4 £) 


— kn 


iMzo) = E (- 1 )^' 


2=0 


{kn + 1 / 2 ) 2 ^ 

j! {2z^oy 


(77) 


In the last relation, /7 means the number of the terms of the asymptotic expansions convenient to take into account. 
We apply asymptotic expansion for the function 2?_fc^_j_i/2 (^o) in Eg. dTSl) . In truncated sum we interchange the 
order of the finite summations over indices k^-i and j. We replace 2? by D, therefore a corresponding power of the 
variable 0 play important role. We obtained the relation: 


^cut 


r(V2) ^ (-1)^ exp(zV4) _kM-- 2 +i /2 

\/27r(l + bA‘^/c)uJo j! {‘^z.'^y ^2 tt( 1 + bA/c) 


E 


fciv-i=0 


V 4uJo / 


(fcAT-l)! 


(78) 


r(fcAr_i + A:jv _2 + 1/2) (/caT-I + 1/2)2^ D-kfi-i-kN- 2 -i /2 {z) 
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Here we simplified the calculations by identity: 

(2fc)! = 22^= fc! (1/2)^, r(fc + 1/2) = r(l/2)(l/2)^ 

Let us study in details the sum 


E V 4i^o > 

{kN-i)\ 


T{kpf-i + kM-2 + 1/2) (JiM-i + 1 / 2 ) 2 ^ ^ —/CAT-l—fcAr-2 —1/2 (■^) 


(79) 


feN'_l=0 


This sum is uniformly convergent, as we could prove by the same procedure as in section ’’Summation over k". 
Therefore we can extend the summation up to inhnity with the error of the order N~^~^. To be able provide the sum 
over index fcjv-i, we must modify term 

(fcjv-i + 1/2)2^ • 

We see, that this object is a polynomial in variable kN-i of the 2j — th order. We rewrite the polynomial in the 
another form: 

min {2j,kN-i) 


{kN-1 + 1 / 2 ) 2 ^ — X/ 


2j (fciv-i)! 


The coefficients can be calculated by recurrence procedure from the relation: 

No 


E 

kN-i—O 


(fejV-l + l/ 2 ) 2 j 
(/cat-i)! 


2i No 

f{kN-i) = ai Y 


1 


2=0 




(/CAT-I - f)! 


TTT f^kn-l) 


From the above definition, we hnd the recurrence equation: 


af = {k-\l2 + i)aY +«^"^ 


(80) 


(81) 


when the initial conditions are: 


= 1 , = (1/2). , aj+i = 0 


The solution of this recurrence equation is: 


a/ = 


( 1 / 2 ),. 

* ) (1/2), 

Inserting all these replacements into Eq. (I75D, with help of the identity 

r(fc/v —2 H“ ^N—\ “h f/2) = r(A;7\r—2 H“ H” 1/2) (/cw —2 i 1/2). 

— 1 — 

after some algebra, by introducing the new summation index k = kjq-i — i we obtain the formula: 


(82) 


^CUt _ 


^/CiV-2 + 1/2 


J 

E 


(-1V 


V4 


V^2(l + 6A2/c)wo ^ /! ( 2 ^ 0 )^ v'27r(l + bA'^/c) ^ 


2j 

E 


(83) 


X r(/c7v—2 H“ ^ H“ 1/2) 


4a;n 


2 \‘^ ^0 ( ^ ^ \ k 
V 4uJn ' 


E 

fc =0 


k\ 


{kN-2 + i + l/2)fe D -few-2—i—1/2-/C 


In the above relation we extend the summation over the index k up to infinity, because the added reminder disappear 
in the continuum limit. The sum over k is now prepared for application of the identity: 


00 / \ 

e"'/^E ^ D-,-k{x) = (x - t) 


k—Q 


(84) 
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The result of the first recurrence step, replacing D by P, can be read: 

1 ^ (-Ip ^2, / 


'p'cut _ __ \'~ i-j _ ^ ^ V -*-/ ^ ^ 

' " V2(l + 5A2/c)u;o ^2^{l + bA^/c)u, ^ j! ( 220 ")^' ^' 


4:U!ooji 


r(fcAr_2 + * + l/2)I?_fc^_2_i_i/2 (^l) 

(85) 


where 

zi = z (1 - ^^/(4wo)) , wi = ^ = 1 - C^/(4wo) ■ 

For the second recurrence step, we have the summation over the index kN -2 in the form: 

1 (^)2fciV-2 


Z 2 - 

kN-2—0 


y/27r(r+"5A27c) (2A:Ar_2)! 


r(A:jv-3 + kN-2 + 1/2)1? — fciV-3 —^JV-2 — 1/2 ('^) 


( 86 ) 


Before replacing in the last equation, let us define the first step of the new recurrence relation by: 


(l)f = 


1 


2i 


LU, 


2i' 


Then the above equation, taking into account the identity zq = zojq in detail can be read: 

1 1 


22 = 


a/ 2(1 + 6A2/c)wo 1 / 2(1 + 5A2/c)wi 

kN-2 


(87) 


X E 

kN-2—0 


1 (^) 

V27r(l + 6A2/c) (2fc^,_2)! + l/2)I?_fc,_3_fc„_,_i/2 (z) 


V ^ V (1)'^' 

(Oy2\j 


j! (2z2)J 

j=0 ' i=0 






r(fcAr_2 + I + l/2)I?_fc„_2_j_i/2 (zi) 


Now, due to the uniform convergence of the sum over index kN -2 we will work with truncated sum of the com¬ 
plete Z 2 . In the truncated sum, we use the asymptotic Poincare type expansion of the parabolic cylinder func¬ 
tion ( 21 ). We have then in the relation for the finite summations only, so we rearrange them, 

in order that summation over kN -2 will be provided first. Remember, that in an intermediate step we replace 
^-fe]v- 3 -fe]v- 2 -i /2 i^) by 7?_fc„_3_fc„_2_i/2 {z) and we must take into account the corresponding power of variable z. 
We have: 


^CUt _ 


1 


1 


1 


a/ 2(1 -|- 5A2/c)a'o a/2(1 -|- bA"^/c)iOi + bA'^/c) 


( 88 ) 


E JziE y (i)2j 

Z^ P (Oy2\j Z^ h 


?! (2z2)t 

i=0 ^ ' i=0 


\ 


UJoOJi 


2 ? = 0 


( 0'' ^ pA/4 kN.2 + 112 


(2z2)' cof 


No 

X E 

kN-2—O 


V4i^l J 


kN-2'- 


-r(fcAr_3 -I- kN-2 + 1 / 2 ) I?_fc„_3_fc;y_2_l/2 (z) (fcjv_2+i/2) 


2l+i 


if the identity 


Z= ^ 
2 


was introduced. Summing over index kN -2 as in the first recurrence step, we have: 

1 1 (W2) 


^cut _ 


y(1 + 6A2/c)u;o y(l + 6A2/c)o;i ^(1 + b^^/c)uj 2 ^ ^ j!n(2z2)^+' a;? 


-fciV-3 .z/_ A/ 

srZr 

j=0 ?=0 


(-l)A' 1 


-27 (89) 


2i 


Ed) 


2i 


i=0 


\ 


i 




LOqUJi 


] E “E* ( - ) r(fcAr-3 +P + l/2)I?_fe„_3_p_i/2(2:2) , 

/ r. \ZJl(+2 J 
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where one postulate the identities: 


By adjustment in the summations: 
J c 


Z2 = ZL02, W2 = 1-■ 


(-iy+^ (-1)^ ^ 


7!Z!(2z 2V+' 
j=0 1=0 ^ ’ 

and interchange of the order of the summations: 


^\{2z- 

/i=0 ^ ' 3=0 


(•••j, 1 = fi-j---) 


2j 2fj.-2j+i fi 2j 

E E -E E 

i=0 p=0 p=0 i=max [0, p—2/i+2j] 


we find the result of the second recurrence step: 

C^'CUt _ ^ 

^2 — 




kN-3 


\/ 2 {l + bA^/c)u!o \/ 2 (l + bA‘^/c)iOi y^ 2 (l + bA‘^lc)oj 2 
J+C, / -r'.p / 42 \P 

^ E ^!(2z2)m E r(A:jv- 3 +P+l/ 2 )X>-fc^_ 3 -p-i/ 2 ( 22 ) 


(90) 


/ j .=0 ^ ^ ^ p =0 

Where the second recurrence step of the function is defined: 


( 2 )^ = E 




1 


3 J UJl 


2p.-2j 


t , (:£r)‘ 


J—0 1 2—max [0, p— 2 /i+ 2 j] 

We can see, that after the A recurrence steps, the result of the A summations over the indices ki can be read: 

r A 


Z cut _ 

A — 


1 


a / 2(1 + bA^fc)u;^ 


(wa) 


— A;a7_i_a 


(91) 


(92) 


Ji-\ - HJa / 2p 


M2z‘^Y 

tl = 0 ^ A J p^Q 

We have postulated the recurrence relations: 


jp 


A2 \ 


WA-lWA 


ZA = ZUJA, WA = 1 - 


jE 

WA-1 ’ 


j r(A:Af-A-i + P + l/2)I7_fcjv-A-i-p-i/2(2;A) ■ 

1/2 + &AVC 


ujo = 


1 + 6A2/C ’ 


and evaluated the recurrence definition for the function (A)p^: 


wT = A 

1=0 


1 


,2p-21 


21 

E 


A —1 i=max [0, p—2^-\-2j] 

when the recurrence procedure begin from: 


A2 


0 -’A-2WA-l 


(A-l)f 


(93) 


(1)?^ = 

Wg 

After the last recurrence step, for A = A — 1, we are left with the relation of the form (15^ where kN-A-i = 0 and 
the index of the V function is only —p — 1/2. We expand the D_p_if 2 izA[-i) as in all previous recurrence steps and 
we find the relation: 


1 


21 


-■N-l 


'N-1 

n 


1 


-HJ/at 


i^o \/2(l + 6A2/c)a;i 


E 

p—O 


E 

1=0 


UJ 


1 

^E 


N-l 


i=0 




UJN-2^N-1 


(-1)^ 

p\{2z^y 

{N-lf/ (1/2) 


2l-\-i 


X 


(94) 
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Following the definition of the al symbols, we have: 


( 1 / 2 ) 


2l+i 


= a 


0 


Then in the last part of the preceding equation we read: 



UJ 


21 

N-1 


21 

■E 

i=0 




CJ]Y-2^N-1 


(iv-i)f 


^2l+i 




(95) 


Following the calculations done in this Appendix, we conclude, that it is possible to provide the summations in 
the precise formula for the N dimensional integral at least by help of the asymptotic expansions of the parabolic 
cylinder functions. It is possible to provide the continuum limit of the our result and there are no additional terms 
contributing to the result in the continuum limit. The result can be read 


Z 


cut 

N-1 



1 

a/2(1 + hAP- lc)uji 




{NfY 


(96) 


This expression is sufficient for calculation of the continuum Wiener unconditional measure functional integral by 
Gelfand-Yaglom procedure leading to the differential equation of the second order. The second part of relation (IMl) 
represent the expansion of an unknown function in present time. 


APPENDIX 2 


In this appendix we will present the idea of the decomposition of the principal result by the summations over 
the indices ki by slightly different method as was presented in the Appendix 1. Our goal is to study the unknown 
function from Appendix 1 in the power expansion in the coupling constant. We start as in Appendix 1 from the 
precise relations for the N dimensional functional integral, and we proceed by the recurrence procedure withouth the 
introducing the function (A)^-^. We start, as in Appendix 1, by summation over the index kM-i- 


^ 1 = E 

fcAr-i=0 


^2 \kN-l 


( — 

\z Zq 


Y^27r(l + bA/c) (2A:Ar_i)! 


r(fc7V-l + Yo) 


-T{kN-2 + ^Af-l + l/2)T’-fc„_2-fe;v-i-l/2 Y) 


Zoy/^ 

Repeating the same calculations as before, we will have after the second recurrence step the relations: 

1 I 


(97) 


^CUt _ 




v'2(l + 5A2/c)a;o v'2(l + bAyc)u;^ ^ 72(1 + 6AVc) iYo 


232 


(98) 






2l=0 


ii 2j2+n 


E 

12—0 


a. 


2j2+n 

'^2 


( 4 . 2122 ) + A + I/2)I?_/c^_3_i2-i/2(2:2) 


In contrary to the method in the Appendix 1, we don’t combine the summations over the asymptotic expansion 
indices. We see, that by this procedure we replace the summations over indices ki by summations over the asymptotic 
expansions indices jA Taking thoroughly into account the last summation term, we find: 


N 


U/2 




Zr = lY[[2{l + bAyc)u;,]\ 

U=1 ) 3l,---,3N 

MM" 


(-ly 


■ N 


E 

2l=0 


*1 J (1/2) 


2ii 


(—) 
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, ( 2 * + V-I ^ ( A- 


X 


2u=0 


2jN-l+lN-2 


-I —2 / Ct • 

E l 2j7V-l + ^N-2 
V iN-i 


(1/2),. 

( 1 / 2 ) 




(99) 


2iV-l=0 


2jiV-l+^iV-2 


( 1 / 2 ) 


22Ar_i 


yl2 


i^7V-2‘^Ar-l 


x(l/2k:„ ,(l/2 + */v-i) 


In the last equation the identities were used: 


where 




4zi-iZi Wi-iWi 




24= = 


and 


By the useful identity 


a? = 


2z 2(l + 6A2/c) 
k\ (1/2), 


* J ( 1 / 2 ), 


= (l/2),^_, (1/2 + V-i),,,. 

we convert our result to the more reliable form for the consecutive calculations: 


N 


- 1/2 


= {l[[2il + bAycU] 

2ji 


L2=1 

2ii 


Jia-aJn 

E 

Jir-'j'iv 

n 


(-ip 


■ AT 


( 1 / 2 ) 


2/1 


ii—O 


E 

2„=0 


2/w-i-|-i«'-2 


2j/i + V—1 


( 1/2 + 


242 


2A+1 


-l-l-iAf-2 / \ / /l2 \ ^ 


*A^-1=0 

In the calculation as follows the key role plays the objects uji defined by recurrence relation 

242 


UJi = 1 - 


(2/i-l 


with the first term 


2/iv 


( 100 ) 


( 101 ) 


U^o — 1/2 + -S 


where 

^ 6A2/C 

2(l + 5A2/c) 

The LOi such defined are represented by the continued fraction. The continued fraction can be represented by the 
simpler relation as the solution of the n — th convergent problem of the continued fraction [H- Let us shortly explain 
this procedure. 
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Let us have a continued fraction of the form: 


Wra — Ol + 


bi 


a2 + 


b2 

03 + - 


The n — th convergent is defined as 


where Pn and qn are defined by equations: 


Qn 


Pn — bln Pn—1 bn Pn—2 
Qn — dn Qn — l T bn Qn—2 


The solutions of these recurrence equations have the form: 


Pn = WiPi + W 2 P 2 
Qn = Wip^ + W 2 P 2 


where pi ^2 are the solutions of the characteristic equation, in our case homogenous one: 

p^ - QnP -bn = 0 

For the continued fraction we have: 

dr). — 1 


( 102 ) 


(103) 


bn = -A^ 

and the solution of the characteristic equation would be: 

Pi ,2 = ^(1 ± V1 - 4^2 ) 

The constants wi and wi are fixed by wq and oji terms, that adjust the initial conditions: 

Po — 1 T 

Pi = l + 2B-A^ 

<70 = 2 

<71 = 1 + 2B 

The n — th convergent method solution is completed by the relations: 

1 


Wl,2 = 


{l + 2B)± v^l - 4^12 + 


2B 


y/l - 4A2 J _ 


W 12 = 1± 


2B 


y/l - 4A^ 

The very important characteristic follows from the above solution: 

Pn — Qn-\-l 

which simplify our calculation significantly. 

We are seeking for an expansion of the functional determinant. It look reasonable to arrange the relation (IIOID 
following the condition 


ji = K = constant 
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where by the value of this constant is labeled the term composed by the sum of all the terms of (IIOII) suitable satisfying 
to the above condition. We can prove immediately, that if this constant, let us choose K, is if = 0 the contributing 
term is only one, it is the relation where all the ji = 0 and value of this contribution is 1. Let us study the nontrivial 
case, when K — 1. The contributions will be all the members of the multiple sum (IIOII) where only one index ji is 
nonzero and equal to one. By some algebra, we can express all the contributions as the sum: 


N 


- 1/2 


P'1 


= \Y[[2{l + bAyc)i0,]\ Sn 


.2 = 1 


(104) 


where 


(ji + i2 + • ■ • + JAT = 1) = 1 + 


-1 


N-l N-1 N-1 


‘^k-l 


/c=l 


'i—k l—n 


y.—k 


1 ! 

m-1 / .2 \ 2 l-l 

\ 


(105) 




l3—m 






In this relation, the term (2 — 5m,i) is due to the binomial coefficients in Ea. IllOlD . Following the result for and the 
functions pi and qi we can calculate the product: 


m —1 

n 

OL—k 




m—1 

n 


A^ 


Qk — lQk 


Pa-l Pa 
— k \ Qa — l Qcjc 


where we have introduced the more convenient variables: 

+”'npj 

and also we performed the replacement: 


Qm—lQr, 


_Ql-i 
Ql 


Inserting the last results into Eg. (11051) we find: 


SnUi + j2 + -h JAT = 1) = 1 - 


( 1 / 2)2 

2^202 


^fc-1 


N-1 

E 

k^l 


• N-1 N-1 

2 E E 

m—k l—m 


Qt-i 


N-1 


Qt-i 

Qm — lQmQl — lQl Q l — lQ I 


-E 

l^k 


In the above equation we have the object symmetric in the indices m and I 

1 


Sm,l — 


Qm — lQmQl — lQl 

This symmetry significantly simplifies the calculation of the double summation, because we can proceed as: 


N-1 N-1 


'N-1 N-1 


N-1 N-1 


Y Sm,l = 1/2 E E E E 


\^m—k l—m 


l—k m—l 


m—k l—m 

In the second term we interchange the order of the summations by identity: 

V —1 N — 1 N — 1 m 


E E = E E 


l—k m—l 


1—k l—k 


Finally, we find: 


N-1 N-1 


'N-1 N-1 


N-1 


E E Sl,m = 1/2 E E Srn,l-Y 


a—k l—r, 


'l—k l—k 


l—k 


(106) 


(107) 


(108) 
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where the identity can be used: 


N-l N-1 


/N-1 

m—k l—k 


. l^k 


Qi-iQi 


To calculate the sum, we introduce the relation: 


Let us remind, that 


where 


Using the identity: 


we are able to calculate the sum: 


Qi = wix* - W2y* . 
Qi = wix* + W2y'' 


x,y = 


Pi,2 


1 


± 


1 

U42 


- 1 


Qi Qi-i 2wiW2{x-y) 


Qi Qi-i 


N-1 

E 

l=k 


Qi—iQi 

( Qn- 1 Qk-1 


Qi-iQi 2 wiW2{x - y) \Qn-i Qk-i 


Inserting these intermediate results into Eg. (11071) . we find: 


(ji + J2 + • • ■ + jw = 1) = 1 — 


( 1 / 2 ), 

2 A2z2 


N-1 


E Qt-i 


k^l 


/N-l 

E 


N-1 


- 2 E 


^ Qi-iQi j Qi-iQi 


The sum over k is the problem of the sum of the finite power series. For the finite N we found the result: 


( ji + ^2 + • • • + jAf — 1) — 1 ~ 


1 


( 1 / 2)2 _ 

2yl2z2 (x — y)^ 


+ 2 


\ x 2N-2 _ y2N-2 

1 1 

[ (x + y) 

x2(x^-l + myN-1^2 y2(yA/-l ma;JV-l)2 

J(l-X2^-2)(1- 



'N-1 I m.yN-l\2 

' W\ ' 


xix^"-^ + ^y 


+ ^x^-i )2 


+ 


(x-y){N -1) 


(x^ ^ 

o .2_ ^ ^ _ 1 

(x^-i + ^y^-^y 


(109) 


( 110 ) 


In the above relation we have all the symbols defined in the previous text. The continuum limit means to take the 
limit IV —> 00 in all formulas obtained for finite N. Doing this, we find for symbols appeared in (IllOl) : 




W2 
Wl 

X — 2 / —>■ 2 A 7 


x,y 

x^ 

-N 


1 


g /37 

C-/37 


A 


2 aA^/c^ 

L 

N 
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Inserting into the Eg. (1110(1 we find finally for the first nontrivial term in the continuum limit the result: 

( 1 / 2)2 « 


SnUi + j2 + ■ • • + JAT = 1) = 1 — 




e2/37 _ e-2/37 


■Pi 


o (e^7_e-d7)2 ■ 

( e ^ T '+ e -^'>')2 


( 111 ) 


In derivation of the above relation no connection between the variables a and b has been supposed. The relation is 
finite in the whole range of the values of 5, positive, zero or negative. For large b and moderate a our result corresponds 
to the results of the conventional perturbative method of calculation. 


APPENDIX 3: THE EVALUATION OF THE FUNCTIONAL INTEGRAL BY THE GELFAND-YAGLOM 

PROCEDURE 

The well-known method of the evaluation of the functional integral for the harmonic oscillator in the continuum 
limit was proposed by Gelfand-Yaglom . We apply the same idea for the evaluation of the functional integral with 
the fourth order term in the potential. 

In the N—th approximation for such functional integral we find in the Appendix 1 the result (1961) and in the Appendix 
2 the approximation of this result (11041) . (11091) . The common form of the both relations for N — th approximation can 
be expressed as: 


2n = 


SN-li^) 


N-l 

I n 2(1 + bA'^/c)uji 

i=0 


The value of the functional integral in the continuum limit is defined as 


Z = lim Zn 

N—*OC 


Let us define the function 


Fk 


n 2(1 -I- bA‘^/c)uji 


i=0 


Sli^) 


( 112 ) 


The function Sk{A) is the result after first k steps od the calculation of the function Sn{A), as it was demonstrated 
in the Appendices 1 and 2. Let us stress that 


Zjq — 


1 


\/Fn 


The aim of the Gelfand-Yaglom construction is to find in the continuum limit the differential equation, such that its 
solution is connected to the continuum functional integral by: 


VW) 

In the sense of the Gelfand-Yaglom construction we are going to express the relation for Fk+i by help of Fk and 
Fk-i- We have: 


Fk+i 


k 

2(1 -I- bA'^/c) n 2(1 -I- bA‘^fc)uJi 


*5^+1 (A) 


k — 1 

n 2(1 -I- bA‘^/c)uJi 


i=0 




(113) 


To obtain the above equation we used (11121) and the identities: 


1 


tUi — 1 — uji_i , ujf) — 1 /2 -)- bZA j c, A 


2(l-h5A2/c) . 
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To introduce Fk and Fk-i, we must to express by help of S'^(A) and S'^_j(A) in the corresponding terms. 


This means, that we use: 


5fe+i(A) = ,5fc(A) + (5fc+i(A)-5fc(A)) 
5fe+i(A) = ,5fc_i(A) + (5fc+i(A)-5fc_i(A)) 


(114) 


In the case of the calculation as in the Appendix 2 the functions S'fe(A) can be calculated without approximations. 
The explicit calculations shown, that (S'/c+i(A) — 5'fc(A))/S'fe(A) A also for the non-perturbative anzatz of the 
Appendix 1. Therefore we expands the denominators in (I113|) up to the terms of the second power in A and we 
obtain: 


1 


1 


5fc+l(A)2 

1 


5fe(A)2 

1 


^ ,5fe+i(A)-5fe(A) , ^ 


/5fc+i(A)-^fe(A)\ 

V Su{A) ) 


+ ... 


5fc+i(A)2 5fe_i(A)2 




^fc(A) 

V ^fc-i(A) 


(115) 


After some algebra we we find for the Eq. (11131) the difference equation: 


Fk+i — 2Fk + Fk-i .Fk — Fk-i Sk+i — Sk 

+ 4--- = i'k 


A 2 


In the continuum limit A 


A A5fc 

0 we use the notation 


2b/i 


_ _ ^ ■S'fc+i — 2Sk + Sk-i _2 [ ‘^fc+i ~ Fk 


A^Sk 


V ASk 


(116) 


kA T 


and we finally find the differential equation: 




d 


d 


_F{r) + 4-E(r) - In S{t) = F{t) S[r) - 4(- In S{r)) 


9r2 


dr 


dr 


2b 




d 


'dr 


(117) 


The initial conditions are: 


m = 1 , 

1 ^(0) = 0- 


By substitution 


F{t) 


y{T) 

S^{r) 


the differential equation can be read in more convenient form: 
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^J/(r) = - y{r) 


(118) 


(119) 


The importance of the thorough evaluation of the function S{t) emerged in this construction of generalized Gelfand- 
Yaglom equation. 


APPENDIX 4: CONDITIONAL WIENER MEASURE 

The conditional measure functional integrals define the propagators in the theory, contrary to the unconditional 
measure functional integrals, which define the partition function. The conditional Wiener measure functional integral 
is defined as the continuous limit of the finite dimensional integral with fixed endpoints: 


1 = 


lim 

A^—»oo 




N 


EA 


c /2 



- <fi-i 

A 



( 120 ) 
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The only difference between the unconditional and the conditional definition inhere in the dimension of the finite 
dimensional integral, whereas the actions are the same. For conditional case, both the endpoints of tp in the time 
variable are fixed. This correspond to calculation of the propagator of the model, whereas the unconditional measure 
integral correspond to calculations of the partition functions. The result of the finite dimensional integration in 
Eg. (11201) . up to index kN -2 is the same as for the conditional measure case. For the last time slice integration, for 
kN-i we have: 


E 

feiv-i=0 


(2A:jv-i)! 


OO 


— OO 


dipN-i 





exp-{aA(p^_i + c/A(l + 5A^/c)<p^_J 


( 121 ) 


The conditional measure result we derive from the result for the unconditional case, when the term depending on the 
summation index k]\[-i we replace by the term: 


cVn 

A(l+6A/c) 


fejv-i 


kM-i=0 


(2fcAr-i)! 


r(fcAr _2 + kpf-i + l/ 2 )I?_fe^_ 2 _f,j^_j_i /2 (z) e 


— {aAc^^ + 




( 122 ) 


In this appendix we do not follow by discussion for the arbitrary value of the ipN, we are going to discuss the 
simplest case of the ip]\f = 0, corresponding to the periodic boundary conditions in the imaginary time variable. The 
result of the N — 1 dimensional integration can be read: 


N-l 


E n 

, • • • — 1 ^ — 1 


iUz) 


2ki 


-\/27r(l + bA/c) (2fci)! 


-F(fci_i + ki + l/ 2 )I?_fe._j_fe^_i /2 (z) 


o~{a^‘PN + 3'(1/2H— 


(123) 


In the last equation the identity fco = 0 is required. The main difference from the unconditional measure calculation 
consists in the fact, that the arguments of the parabolic cylinder functions in the Ea. (ll23() are the same in contrary 
to unconditional measure case, where the argument of the last function differs from the all others arguments of the V 
functions. We are going to repeat the calculation as in the Appendix 2, to show the differences of the unconditional 
measure calculations from conditional measure ones for the first two terms of the decomposition of the corresponding 
functional integrals. We again introduce the objects as for unconditional measure case, but corresponding to the 
conditional measure specifications: 


A2 

UJi = I - 

UJi-1 

with 

Wo = 1 . 

Following the unconditional measure calculation, by the n-th convergent method we find: If w„ is the continued 
fraction defined by: 




Cll + 


hi 


a2 + 


asH— 


then we can find the ojn in the form of the n — th convergent: 




Qn 


where pn and Qn are the solutions of the homogeneous difference equations with the constant coefficients: 


Pn — Pn—1 H" Pn—2 
Qn — Qn — 1 Qn—2 

(124) 
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The solutions have the form: 


Pn = WiPi + W 2 P 2 


q-n = WlPi + W 2 P 2 

where pi ^2 are the solutions of the characteristic equation: 

- QnP -b„ = 0 

where 

— 1 

bn = -A\ 

The boundary conditions for are fixed by the first two relations for u;„: 

Wo = 1 

wi = 1 - 

and we have: 

Po = 1 
Pi = 1-A^ 

qo = 1 
qi = 1 

The n — th convergent solution for w„ is thus given by the relations: 

1 


(125) 


Pi ,2 = -(1 ± a/i - 4^2 ) 

1 / 1-2^2 
W1.2 = - I 1 ± 


(126) 

(127) 

(128) 

(129) 

(130) 

(131) 

2 VI-4.42; 

We see from the last relations the very powerful characteristic of the solution represented by the identity: 

Pn = qn+i, n = 1, 2, • • • 

Expanding the — 1 dimensional integral for the conditional measure into the asymptotic series, we find the 
relation: 

(N-l 'j ~^/2 

= <^ n [2(1 + ^AVc)w,] \ Zn (133) 


Vl - 4^12 

1 / 1 

Wl,2 = - 1± 


i=l 


The first part of this relation in the continuum limit is the conditional measure functional integral determinant for 
the harmonic oscillator type action. The second part, marked Zpj, is the asymptotic expansion correction. As for the 
unconditional measure integral we calculated the first two corrections terms and we find: 

(^2 


1 - 


1 


2z'^ {x — y)^ 


(134) 


[ a.2JV-2 _ y2N-2 

1 1 

1 

\ {x + y) 



+ 2 


(1 I 1£2_ _ rr'2'N-2 _ W2,,2N-2\ n _i_ Hi _ i,2A''-2 _ wXrr‘^N-2\ 

V roi bUi ^ ^ ^ Wo y Wo / 




y{y^-^ + ^x^-i)2 


+ 


{x-y){N - 1) 


/ iV-1 _ W2 N-l\2 
o ^ ^ _ 1 
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Above relation possesses the well defined continuum limit. We find, in the limit N —>■ oo the following relations: 

-1 


W2 

Wi 

x-y 




2A7 

1 

g/37 

p-01 


1 

-2 

A 


2aA^/c^ 

>L 

N 


then, finally, in the continuum limit the first nontrivial correction to Zn is: 

1 


1 - 


(5)2“ 


873 - e-'yf^y 


{_3(e27/3 _ e-27/3) + 27/3 + 4] } 


(135) 


or, 


1 - 


3a 

320^73 


3 coth(7/3) + 27/3 


coth^(7/3) + 


1 


2sinh (7/3)J 


(136) 


We can stress, that our result is valid for finite values of the parameter ”6” positive or negative, so we describe the 
case of the anharmonic oscillator and the Higgs model as well. 
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